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Abstract- In radial distribution systems, the voltages at buses 
reduces when moved away from the substation, also the losses 
are high. The reason for decrease in voltage and high losses is 
the insufficient amount of reactive power, which can be 
provided by the shunt capacitors. For this purpose, in this 
paper, two stage methodologies are used. In first stage, the 
load flow of pre-compensated distribution system is carried 
out using 'Dimension reducing distribution load flow 
algorithm'. In the second stage, Differential Evolution (DE) 
technique is used to determine the optimal location and size 
of the capacitors. The above method is tested on IEEE 69 bus 
system. In this paper a new method is proposed to improve the 
power factor of those buses having low power factor (less than 
0.81ag) to unity power factor simultaneously by placing the 
capacitors. 

Index Terms- Electrical Distribution System, Optimal 
Capacitor Placement, Dimension reducing distribution load 
flow algorithm (DRDLFA), Differential Evolution (DE) 
Algorithm, power factor correction. 

I. Introduction 

Capacitors are generally used for reactive power 
compensation in distribution systems. The purpose of 
capacitors is to minimize the power and energy losses and to 
maintain better voltage regulation for load buses and to 
improve system security. The amount of compensation 
provided with the capacitors that are placed in the distribution 
network depends upon the location, size and type of the 
capacitors placed in the system [1]. A lot of research has 
been made on the location of capacitors in the recent past [2] 
without including the installation cost of the capacitors. All 
the approaches differ from each other by the way of their 
problem formulation and the problem solution method 
employed. Some of the early works could not take into 
account of capacitor cost. In some approaches the objective 
function considered was for control of voltage. In some of 
the techniques, only fixed capacitors are adopted and load 
changes which are very vital in capacitor location was not 
considered. Other techniques have considered load changes 
only in three different levels. A few proposals were schemes 
for determining the optimal design and control of switched 
capacitors with non-simultaneous switching [3]. It is also 
also very important to consider the problem solution methods 
employed to solve the capacitor placement problem, such as 
gradient search optimization, local variation method, 
optimization of equal area criteria method for fixed capacitors 
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and dynamic programs [4]. Although these techniques have 
solved the problem, most of the early works used analytical 
methods with some kind of heuristics. In doing so, theproblem 
formulation was oversimplified with certain assumptions, 
which was lacking generality. There is also a problem of local 
minimal in some of these methods. Furthermore, since the 
capacitor banks are non continuous variables, taking them 
as continuous compensation, by some authors, can cause 
very high inaccuracy with the obtained results. A differential 
evolution algorithm (DEA) is an evolutionary computation 
method that was originally introduced by Storn and Price in 
1995 [5] . Furthermore, they developed DEA to be a reliable 
and versatile function optimizer that is also readily applicable 
to a wide range of optimization problems [6] . DEA uses rather 
greedy selection and less stochastic approach to solve 
optimization problems than other classical EAs. There are 
also a number of significant advantages when using DEA, 
which were summarized by Price in [7]. Most of the initial 
researches were conducted by the differential evolution 
algorithm inventors (Storn and Price) with several papers [8] 
which explained the basis of differential evolution algorithm 
and how the optimization process is carried out. In this 
respect, it is very suitable to solve the capacitor placement or 
location problem. IEEE 69 bus distribution system is 
considered for case study. The test system is a 12.66 KV, 10 
KVA, 69-bus radial distribution feeder consisting of one main 
branch and seven laterals containing different number of load 
buses. Buses 1 to 27 lie on the main branch. Bus #1 represents 
the substation feeding the distribution system. 

II. Distribution Load Flow 

The distribution systems are characterized by their 
prevailing radial nature and high R/X ratio. This renders the 
load flow problem ill conditioned. So many methods [9-14] 
have been developed and tested ranging from sweep 
methods, to conic programming formulation. The relationship 
between the complex branch powers and complex bus powers 
is derived as a non singular square matrix known as element 
the equations are expressed in matrix format. This proposed 
dimension Reducing load flow method could be applied to 
distribution systems having voltage-controlled buses also. 
Notations 

N-no of buses 

I ij -Branch current flowing through element ij 
I j -Bus current of node j 
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V j -Bus voltage of node j 

S ij -Complex power flowing from node i to node j 
S ji-Complex power flowing received at node j from 
node i 

S j -Specified Bus power at bus j 
Z ij -Impedance of element ij 
TL ij-Transmission loss of element ij 
The power flow method is summarized as follows: 

1 . For the first iteration transmission losses are initialized 
as zero for each element. 

2. From the bus powers specified the branch powers are 
determined as per equation (1&2). 



J - 1T~ l T 
1 branch ~ ^ 1 bus 

ci _ z*\o sending _ j-t 



(1) 

(2) 



3. The variable Rij. is determined for each element using 
equation 3 



: P. + l O, 
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(3) 

4. The bus voltage, branch current and bus current are 
determined from Rij. 

K 

V^V-t^ (4) 



v: 
i 

5. The bus currents are determined 



(5) 



(6) 



And from (6) bus powers are calculated. Since the 
transmission losses are neglected in the first iteration there 
will be mismatch between the specified powers and calculated 
powers. The mismatch is a part of the transmission loss. TLijr 
is the transmission loss part for 'ij'th element for 'r'th iteration. 
Transmission loss of each element is the summation of the 
transmission loss portions of all previous iterations. 



TL S .=^\ 'TL „ 

Where is the Iteration count 
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max( TLV} < £ 
Treatment of voltage controlled buses If power is fed from 
multiple ends of the radial system, other feeding buses except 
slack bus are treated as voltage controlled buses. The 
equation is as follows. 



(10) 



The trigonometric equations are to be solved to get the phase 
angle of each PV bus j and the reactive power can be updated 

as 



Equation 10 is modified for the jth voltage controlled bus. 
real(S t ^) = R = ^al(R^) (12) 
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Then the same procedure described for the PQ buses is carried 
out till the convergence. 

III. POWERFACTOR CORRECTION 

By supplying leading KVAR with the addition of the 
capacitor, the total lagging KVAR required by the load can 
effectively be reduced and the apparent power S of the load 
can be reduced. This effectively decreases the power factor 
angle and improves the power factor. The generation of 
reactive power at a power plant and its supply to a load 
located at a far distance is not economically feasible, but it 
can easily be improved by capacitors located at the load 
centers. 

Assume that a load is supplied with a real power P, lagging 
reactive power Q , and apparent power Sj at a lagging power 
factor of 

coe6i = P1/S1; where Si = (P 2 + Cb) 1 ' 2 (Y1) 

When a shunt capacitor of Q is installed at the load, the 
power factor can be improved from cos8, to cos8„ 

where = [P 2 +( Qi-Q c ff 2 ; QiKfc+Cfc (18) 

Thus the power factor correction produces economic savings 
in capital expenditures and fuel expenses through a release 
of kilovolt ampereage capacity. 

a. Savings Calculations Due To Power Factor 
Correction 
•PSERate Schedule: 
• Energy Rate = $4.08 per KWH 
•Energy Rate = $4.08 per KWH 

IV. Differential Evolution 

Differential Evolution is relatively a newer addition to the 
population based search algorithms. DE was first suggested 
by Storn and Price in 1 995 as a search technique for solving 
optimization problems. It uses the same operators like 
mutation, crossover and selection as that of GA but 
manipulates them in a manner different to that of GA. 
Differential Evolution (DE) has proven to be a promising 
candidate for minimizing real valued, multi-nodal objective 
functions. Besides it has good convergence properties. DE 
is very simple to understand and to implement. 
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Differential evolution (DE) is a population-based 
stochastic optimization algorithm for real- valued optimization 
problems. In DE each design variable is represented in the 
chromosome by a real number. The DE algorithm is simple 
and requires only three control parameters: weight factor (F), 
crossover rates (CR), and population size (NP). Three main 
steps of DE, mutation, crossover, and selection were 
performed sequentially and were repeated during the 
optimization cycle. Figure 1 gives the flowchart of DE. 
a. De Implementation 

The basic procedure of DE is summarized as follows. 

Step 1 : Randomly initialize population of individual for 

DE. 

Step 2: Evaluate the objective values of all individuals, 
and determine the best individual. 

Step 3: Perform mutation operation for each individual in 
order to obtain each individual's corresponding mutant 
vector. 

Step 4: Perform crossover operation between each 
individual and its corresponding mutant vector in order to 
obtain each individual's trial vector. 

Step 5: Evaluate the objective values of the trial vectors. 

Step 6: Perform selection operation between each 
individual and its corresponding trial Vector so as to generate 
the new individual for the next generation. 

Step 7: Determine the best individual of the current new 
population with the best Objective value then updates best 
individual and its objective value. 

Step 8: If a stopping criterion is met, then output gives its 
bests and its objective value. 

Otherwise go back to step 3. 



Initial the 

DE 
Population 



Calculate the fitness of all 
population 



Mutation 
Crossover 
Selection 




[ Stop ) 



Fig 1: Differential Evolution flowchart 



V. Mathematical Formulation 

The cost function is formulated and it can be represented 
by the equation 3.1 [15]. 
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> (19) 



Total cost in $ 

Factor to convert energy losses to dollars = 

$0.06/KWH 

Load Duration in Hours 

Real Power loss in KW 

Capacitor Cost = $ 3/KVAR 

Reactive power injection from capacitor to node i 

Capacitor Installation Cost = 1000$/KVAR 

Number of Capacitor locations 

VI. Case Study And Results 



Prior to capacitor installation, a load flow program based 
on Dimension Reducing Distribution load flow method is run 
to obtain the present system conditions. The proposed solu- 
tion methodologies have been implemented in MATLAB 
7. 10.0. The solution algorithm is based on DE algorithm and 
tested on IEEE 69 Bus System in Fig. 2 which has been de- 
signed to find the optimal solution for this problem. 
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Fig 2: IEEE 69 Bus System 

The parameters are defined as shown below: 

Elapsed time: 92.38 Sec; Gmax = 800; F: 0.8; CR: 0.8; NP: 
100 
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Table 1 : Power Factor Correction 



Bus No. # 


Old p.f 


New p.f 


Capacitor rating 
required in KVAR 


5 


0.7633 


1 


2.2 


9 


0.7791 


1 


3.5 


40 


0.768 


1 


1.0 



Table 1 shows that there are three buses having low power 
factor i.e less than 0.8 out of 69 buses and the buses are #5, 
#9 and #40. The rating of the additional capacitors required 
for power factor correction to unity power factor 
simultaneously at buses #5, 9 and 40 are 2.2, 3.5 and 1.0 
KVAR respectively. 



Table II: Savings Due To Power Factor Correction 





Bus No. #5 


Bus No. #9 


Bus N). #40 


Savings inEhergy 
usage inVyear 


193884 


3054.72 


879.12 


Savings duetopf. 
inprovernat in 
Vjear 


2890.8 


4599 


1314 


Tctal savings in 
Vjear 


482964 


7653.72 


2193.12 



Table 2 gives savings that can be obtained by correcting 
the power factor to unity at the buses #5, #9 and #40 by 
placing the additional capacitors. Table 3 gives the ratings 
of the capacitors required only for improving the voltage 
profile and reducing the losses using DE. 



Table III: KVAR required only for improving the voltage profile and 

REDUCING LOSSES USING DE 



Optimal Location 


Optimal Size of the capacitor 


Bus No #16 


200 KVAR 


Bus No #60 


700 KVAR 


Bus No #61 


500 KVAR 



Table 4 gives the total KVAR capacity required including 
additional capacitors to improve the power factors of the 
weak buses to unity in addition, to improve the voltage 
profile and to reduce the energy losses using DE on IEEE 69 
Bus System. 

Table IV: Total KVAR placed including additional capacitors on IEEE 
69 Bus System for p.f.c using DE 



Optimal Locati en 


Optimal Sizeof the capacitor 


Bus No #16 


200KVAR 


Bus No #60 


700 KVAR 


Bus No #61 


500 KVAR 


Bus No #5 


2.2 KVAR 


Bus No #9 


3.5 KVAR 


Bus No # 40 


1.0 KVAR 



Table 5: Comparison of Per Unit voltages of IEEE 69 Bus System of 
DE method before and after p.f. correction 



Bus No. # 


Voltage in (p.u.) 
before placing 
capacitors 


Voltage in 
(p.u.)uaneDE 
afterplacing 
capacitor 


Voltage in 
(p.u.) using DE 
after p.f 
conection 


1 


1.0000 


1.0000 


1.0000 


2 


1.0000 


1.0000 


1.0000 


3 


0.9999 


1.0000 


1.0000 


4 


0.9998 


0.9999 


0.9999 


5 


0.9990 


0.9994 


0.9994 


6 


0.9901 


0.9923 


0.9923 


7 


0.9808 


0.9849 


0.9849 


8 


0.9786 


0.9831 


0.9832 


9 


0.9774 


0.9823 


0.9823 


10 


0.9724 


0.9776 


0.9777 


11 


0.9713 


0.9766 


0.9766 


12 


0.9682 


0.9738 


0.9738 


13 


0.9653 


0.9713 


0.9713 


14 


0.9624 


0.9689 


0.9689 


15 


0.9595 


0.9665 


0.9665 


16 


0.9589 


0.9661 


0.9661 


17 


0.9581 


0.9653 


0.9654 


18 


0.9581 


0.9653 


0.9654 


19 


0.9576 


0.9649 


0.9649 


20 


0.9573 


0.9646 


0.9646 


21 


0.9568 


0.9641 


0.9641 


22 


0.9568 


0.9641 


0.9641 


23 


0.9567 


0.9640 


0.9640 


24 


0.9566 


0.9639 


0.9639 


25 


0.9564 


0.9637 


0.9637 


26 


0.9563 


0.9636 


0.9636 


27 


0.9563 


0.9636 


0.9636 


28 


0.9999 


0.9999 


0.9999 


29 


0.9999 


0.9999 


0.9999 


30 


0.9997 


0.9998 


0.9998 


31 


0.9997 


0.9997 


0.9997 


32 


0.9996 


0.9996 


0.9996 


33 


0.9993 


0.9994 


0.9994 


34 


0.9990 


0.9990 


0.9990 


35 


0.9989 


0.9990 


0.9990 


36 


0.9999 


0.9999 


0.9999 


37 


0.9997 


0.9998 


0.9998 


38 


0.9996 


0.9996 


0.9996 


39 


0.9995 


0.9996 


0.9996 


40 


0.9995 


0.9996 


0.9996 


41 


0.9988 


0.9989 


0.9989 


42 


0.9986 


0.9986 


0.9986 


43 


0.9985 


0.9985 


0.9985 


44 


0.9985 


0.9985 


0.9985 


45 


0.9984 


0.9984 


0.9984 


46 


0.9984 


0.9984 


0.9984 


47 


0.9998 


0.9998 


0.9998 
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Table 5: Comparison of Per Unit voltages of IEEE 69 Bus System of 

DE METHOD BEFORE AND AFTER P.F. CORRECTION 



48 


0.9985 


0.9986 


09986 


49 


0.9947 


0.9948 


09 948 


50 


0.9942 


0.9942 


09 942 


51 


0.9785 


0.9831 


09 831 


52 


0.9785 


0.9831 


09 831 


53 


0.9747 


0.9802 


09 802 


54 


0.9714 


.9 779 


09779 


55 


0.9669 


0.9746 


09 746 


56 


0.9626 


0.9715 


09715 


57 


0.9401 


9 538 


09 538 

{J.y J JO 


58 


0.9290 


0.9451 


09 451 


59 


0.9248 


0.9417 


09417 


60 


0.9197 


.9 377 


09377 


61 


0.9123 


0.9325 


09326 


62 


0.9120 


.9 324 


09324 


63 


0.9117 


0.9321 


09321 


64 


0.9098 


0.9302 


09302 


65 


0.9092 


0.9296 


09296 


66 


0.9713 


0.9766 


09766 


67 


0.9713 


0.9766 


0.9766 


68 


0.9678 


0.9735 


0.9735 


69 


0.9678 


0.9735 


09735 



Table VI: Savings due to power factor correction of IEEE 69 Bus 
System using DE 





Before 
placing 
capacitor 


After placing the 
capacitor* before 
p.f. correction 


After placing 
the capacitor& 
after p.f. 
correction 


Energy loss cost/ 
Year 


1,18,260$ 


77,767.77$ 


77,741.49$ 


Saving due to 
reduction in Energy 
loss / year 




40,492.23$ 


40,518.51$ 


Savings duetopf. 
correction /year 






14,676.48$ 


Total cost of the 
Capacitors (KVAR 
cost) 




4200$ 


4220.1$ 


Installation cost of 
the Capacitors 




3000$ 


6000$ 


Net Savings / Year 




33,292.23$ 


44,974.89$ 



Table 5 gives the comparison of Per Unit voltages of 
IEEE 69 Bus System before and after p.f. correction and 
Table 6 gives the savings due to power factor correction of 
IEEE 69 Bus System using DE. Finally the net savings before 
p.f. correction are 33,292.23$ and 44,974.89$ per annum 
after p.f. correction. 

Table 7 gives comparison between savings due to power 
factor correction of IEEE 69 Bus System using DE and GA 
[15] . The reduction in real power losses is almost same in I. 
in both GA and DE. The total number of capacitors required 
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using GA is 7 whereas for DE it is 6. Because of this the 
reduction in the installation cost of the capacitors is $1000 
and also the total KVAR required in DE is 1406.7KVAR 
whereas in the case of GA it is 1506.7KVAR. Because of this 
also the reduction in the cost of the capacitors is $300. 
Therefore when compared to GA, DE is giving better savings. 
DE is the best method when compared to GA because of its 
savings, less number of capacitors to be installed and 
reduction in the installation cost of the capacitors even 
though there is slight improvement in the magnitudes of the 
voltage profile in GA when compared to DE. 

Table VII: Comparison between savings due to power factor correction 
of IEEE 69 Bus System using DE and GA 





After 
placing 
capacitor & 
before p.f .c 
using 
CA[15] 


After placing 
capacitor & after 
p.fc. 
using CA[ 15] 


After placing 
capacitor& 
before p.f c. 
using DE 


Afterplacing 
capacitor & after 
p.f.c. 
Using DE 


Energy loss 
cost/Year($) 


77,152.82 


77,131 J8$ 


77,767.77 


77,741.49 


Sa\ingsdue 
Energy loss 
reduction/ 
year($) 


41,107.18 


41,1282$ 


40,49223 


40,518.51 


Sa\ingsdue 
topf.c/ 
year($) 




14,676.48$ 




14,676.48 


Total 
Capacitors 
cast® 


4,500 


4,520.1 


4,200 


4,220.1 


Installation 
cast® 


4,000 


7,000 


3,000 


6,000 


Net Savings/ 
Year($) 


32,607.18 


44,284.58$ 


33,292.23 


44,974.89 



vn. CONCLUSION 



Buses having low power factor (less than 0.8 lag) are 
identified for the LEEE 69 Bus System and their power factor 
had been improved to unity simultaneously with the help of 
additional capacitors by DE Algorithm along with the 
Dimension Reducing Load Flow method. The act of improving 
the power factor of the poor power factor buses to unity 
power factor has resulted in better voltage profile, reducing 
losses in spite of additional capacitors placed, due to further 
increase in savings because of power factor correction. In 
overall consideration the Differential Evolution algorithm 
method along with improvement of the poor power factor 
buses to unity power factor was found to give good voltage 
profile, reduced losses and improved savings from the point 
of view of annual savings with good performance this seems 
to be the Best Method. 
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